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From observations of the metallic species sodium (Na), potassium (K), and magnesium (Mg) in Mercury's 
exosphere, we derive implications for source and loss processes. All metallic species observed exhibit a 
distribution and/or line width characteristic of high to extreme temperature - tens of thousands of 
degrees K. The temperatures of refractory species, including magnesium and calcium, indicate that the 
source process for the atoms observed in the tail and near-planet exosphere ale consistent with ion sput¬ 
tering and/or impact vaporization of a molecule with subsequent dissociation into the atomic form. The 
extended Mg tail is consistent with a surface abundance of 5-8SS Mg by number, if '30% of impact-vapor¬ 
ized Mg remains as MgO and haif of the impact vapor condenses. Globally, ion sputtering is not a major 
source of Mg, but locally the sputtered source can be larger than the impact vapor source. We conclude 
that the Na and K in Mercury’s exosphere can be derived from a regolith composition similar to that of 
Luna 16 soil (or Apollo 17 orange glass), in which the abundance by number is 0.0027 (0.0028) for Na and 
O.O006 (0.0045) for K. 
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t. Introduction 

We present observations of the metallic species sodium (Na). 
potassium £K), and magnesium (Mg) in Mercury’s exosphere. 
Although sodium and potassium have been observed previously 
from ground-based platforms (e.g.. Potter et al„ 2002), we present 
a map of the Na/K abundance ratio over the observed sunlit disk of 
the planet. Magnesium was discovered in Mercury’s exosphere 
during the second flyby of the innermost planet by the MErcury 
Surface, Space ENvironment. GEochemistry, and Ranging (MES¬ 
SENGER) spacecraft (McClintock et al., 2009). Here, we present 
an analysis of the observed Mg column abundances and a discus¬ 
sion of the implications for source and loss processes. All metallic 
species observed exhibit a distribution and/or line width charac¬ 
teristic of high to extreme temperature. The temperatures of 
refractory species, including magnesium (Mg) and calcium (Ca), 
indicate that the source process for the atoms observed in the tail 
and near-planet exosphere are consistent with ton sputtering and/ 
or impact vaporization of a molecule wi th subsequent dissociation 
into the atomic form. The very high energy metaiiic species may 
indicate that a resonance process of some kind is acting to produce 
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a highly energetic gas that may not be thermal. We discuss the re¬ 
sults in the context of assumed surface composition, and we assess 
the source and loss processes as constrained by the modeled mag- 
netospheric configuration at the time of observation and other 
ancillary information. 


2. Observations 

We report in Section 2.1 observations of potassium and sodium 
in Mercury's exosphere taken at the McMath-Pierce solar tele¬ 
scope on Kitt Peak, Arizona, in support of the first MESSENGER fly¬ 
by (Ml) on January 14. 2008. In Section 2.2, we present 
observations of magnesium taken with the Ultraviolet and Visible 
Spectrometer (UWS) channel of MESSENGER’S Mercury Atmo¬ 
spheric and Surface Composition Spectrometer (MASCS) instru¬ 
ment during the spacecraft's second Mercury flyby (M2) on 
October 6, 2008 (McClintock et al„ 2009). Although potassium 
has been observed previously (Potter and Morgan. 1997: Potter 
et a!., 2002, and references therein), we present a map of the ratio 
of sodium to potassium in Mercury’s exosphere. Physical data per¬ 
tinent to the observations or data reduction are listed in Table t. 
The high Doppler shift of Mercury with respect to the Sun is advan¬ 
tageous for these observations. 
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Table 1 

Mercury ephemerides for observations. 

Date Soec.es H rA (AU) r, heliocentric Phase Angular 

radial velocity aI2 g!e {'■) diameter 
_{kms-'j_(ataecj 


January 17.2006 Na, K 0336 -8.72 63 6.2 

October 6,2008 Mg 0.342 -9.16 173 Hi A 


2.1. Observations of the sodium and potassium exospheres 

Observations of the potassium exosphere of Mercury were ob¬ 
tained at the McMath-Pierce solar telescope on January 17, 2008, 
three days after the first MESSENGER flyby. We observed the K 
D 2 line at 766,49 nm at rest in air; the K D, line is at 
769.8974 nm in air. A large Doppler shift is advantageous for 
potassium observations in order that the atomic line at rest with 
respect to the planet should coincide with the continuum region 
of the solar flux and be removed from the telluric absorption lines 
as illustrated in Fig. 1. We used a 5" x 5" image slicer (where " de¬ 
notes arcsec) to direct the light from a 5” x 5" square area on the 
planet to the entrance slit of the stellar spectrograph at the tele¬ 
scope. Because the angular diameter of Mercury as viewed from 
Earth was 6.17", three different observations covering a 5" x 5" 


area were mosaicked together to produce a 5" x 10" image cover¬ 
ing the entire illuminated planet (Fig, 2). The phase angle of Mer¬ 
cury as seen from Earth was 63°, and we were primarily 
observing the dawn side. 

The K observations were reduced to column abundance using a 
K D-i g-value (photons s 1 atom 1 ) of 107.7 computed at a helio¬ 
centric radial distance R orb »0.336 AU and a heliocentric radial 
velocity, R ~ 8.7 km s 1 (16:00 UTC). Note that the solar flux as 
seen at Mercury is shifted blueward (toward shorter wavelengths), 
but the emission in the exosphere is on the fed side of the solar 
Fraunhofer feature as seen from Earth. The g-vaiue (pho¬ 
tons s - ' atom 1 ) was computed using the solar flux atlas at theKitt 
Peak National Observatory (KPNO) McMath-Pierce solar data ar¬ 
chives (Kurucz et al„ 1984) and an oscillator strength./-value, of 
0.682 for K D 2 from the National Institute of Standards and Tech¬ 
nology (NIST) atomic database (Ralchenko et a!., 2008). In spite 
of the large g-value and high solar continuum, the line is still very 
weak. Pertinent observational data for the K line are listed in 
Table 2. The heliocentric distance of Mercury. Rorb. was obtained 
from the jet Propulsion Laboratory (JPLj Horizons website at the 
mean time of the observations. The photo-ionization lifetime of K 
listed in Table 2 is from Fuiie et al. (2007); that for Na is derived 
from the theoretical cross section (Huebner et a!., 1992), both at 
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Table 2 

Physical data for the potassium and sodium observations. 


Wavelength (nm ), K (D,5 (rest) 
g-value K(D 2 ) 

K photo-ionization lifetime 
Wavelength {nm}, Na (Dj) 
g-value Na (D 2 ) 

Na photo-ionization lifetime 
Sub-solar longitude 
Sub-Earth longitude 
Mean anomaly 


Value 


766-49 

107.72 

5.1 x» 3 s{T.4h) 

586.895 

38.2 

1.86 * (5.2 h) 

1S2“£ 

120°E 


fiotb = 0.336 AU. Because the photo-ionization rates depend on the 
EUV flux, which depends on the solar cycle, we have used the rates 
calculated for 'quiet' Sun. 

We obtained observations of the sodium D lines with the 
McMath-Pierce solar telescope using the 5" x 5" image slicer on 
each day during January 12-17, 2008, in support of the first MES¬ 
SENGER flyby. The Na D 2 line at 588.995 nm (at rest) was observed 
with the same instrument and identical setup as for potassium ex¬ 
cept for the wavelength setting. The data reduction process has 
been described in detail by Kiilen etai, (2001), and the calibration 
procedure has been described by Kiilen et al. (1990). Because the 
image of Mercury was slightly larger than the field of view of the 
image slicer, mosaicked images of both sodium and potassium 
were made to determine the map of Na/K emission above the sun¬ 
lit disk of Mercury (Fig. 3). Note that the sodium D line emission is 
much more intense than the potassium D line emission, as illus¬ 
trated in Fig. 1, primarily because the K abundance is much less 
than that of Na, Both the background solar flux and the oscillator 


strength of the potassium D Sines are strong. The signal to noise 
(S/N) in the Na D 2 emission is 100 for the brightest pixel in the fig¬ 
ure. The S/N for the K 0 2 emission tine is 7 for the brightest pixel. 
Thus the noise-equivalent column is 10® cm 2 for Na and 7 x 10* 
for K. All of our observed column abundances in Na had S/N values 
between 100 and 10, and our S/N values for K were between 10 and 
1.4, The K abundance in the northern hemisphere is very low, espe¬ 
cially at high latitudes. 

Although the pixels in the K data shown in Fig. 2 were averaged 
over a 2 x 2 pixel region, the actual observed pixel elements were 
retained so that the division could be done. Potassium pixels in the 
fourth and fifth rows from the bottom (south) and columns 1-8 
from left (west) were replaced by the averages of rows 2, 3, 6. 
and 7 in their respective columns (i.e„ the adjacent rows). Those 
pixels appeared to have data dropouts (see Fig. 2) but were off 
the planet except for the last two and would not contribute to 
the average Na/K ratio. The iine-of-sight Na/K ratio in the southern 
hemisphere (~36) is dose to the smallest ratio found in the data¬ 
sets considered by Potter et al. (2002), whereas the ratio in the 
northern hemisphere (—74} is near the mid-value found previ¬ 
ously. Note that Potter et al, (2002) published only the average of 
the pixel ratios, although the ratios were computed in the same 
way as for Fig. 3A. The large Na/K ratio in the northern hemisphere 
reflects the fact that most of the K values in the northern hemi¬ 
sphere are noise-equivalent column. If we exclude from the aver¬ 
age those pixels for which the K is at noise-equivalent column 
and ratio the zenith column abundances, the Na/K ratio in the 
northern hemisphere is 49 and that in the southern hemisphere 
is 22. There is a small excess of Na in the northern hemisphere rel¬ 
ative to the global average value and a large excess of K in the 
southern hemisphere. The pixels outside of the illuminated disk 



Fig. 1 A (left) Sodium-co-potassjum ratio map from Hne-of-sigitt column abundances. B (right) Ratio from column abundances reduced to zenith columns. The values of Na/K 
ratio greater than 70 are dominated by noise inK, Ratios In the southern hemisphere and those near the equator in the northern hemisphere of-25-35 are obtained with K 
observations having signal/noise >5. The Na/K in the northern hemisphere computed from the Iine-of-sight ratios is 74, and that in the southern hemisphere is 36. The sub- 
Earth point is marked with an x. North is at the top. and the terminator is to the right. 
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were not counted in the ratios because the K line is weak and noise 
dominates the ratio off the disk. In addition to noise, the difference 
in scale heights would affect the ratio off the disk. 

Zenith column abundances, the integrated number density nor¬ 
mal to the surface, were derived for each pixel from the line-of- 
sight column abundances at the same pixel using a Chamberlain 
exosphere and a temperature of 1S00K for both Na and K 
(Chamberlain and Hunten, 1987, chapter 7), under the assumption 
of a locally uniform exosphere. The derived average Na zenith col¬ 
umn in the north, 4.4x10 ,0 cm 2 , and that in the south, 
3.3 x 10'° cm 2 , are within the range previously derived by Potter 
et al. (2002), taking into account that Potter et al. (2002) published 
iine-of-sight abundances. The average did not include the central 
row (roughly equatorial) because there was an odd number of rows 
on the planet disk, but the addition of this row did not change the 
averages. The average K zenith column in the north is 8 x 10 s and 
that in the south is 13 x 10 s . The derived average K abundance in 
the north is very close to the noise-equivalent column (7 x 10 s ), 
while that in the south is 1.9 times the noise-equivalent column. 
The average Na abundance in the south divided by the average K 
abundance in the south is 22, excluding pixels for which K is 
noise-equivalent column, and this is the same as the average of 
the Na/K computed per pixel. The average of the ratios of Na/K 
computed per pixel in the north is dominated by the values near 
the limb where the K column is the noise-equivalent column un¬ 
less those pixels are excluded from the average. The Na spatial var¬ 
iation represents only a 10% excess in the north relative to the 
average value over the entire Earth-facing illuminated disk. The 
observed K that is above noise level is in the southern hemisphere 
or close to the equator in the northern hemisphere. This local var¬ 
iation is different from the global variation seen in the datasets 
analyzed by Potter et al. (2002), who found that K and Na vary sim¬ 
ilarly (i.e„ both north or south). We find that the Na excess is in the 
northern hemisphere and the K excess is in the southern hemi¬ 
sphere. The Na is more or less peaked at the subsolar point. How¬ 
ever, the distribution of K is different. In the northern hemisphere, 
K is almost entirely confined to within 30“ of the equator and near 
the subsolar point, where the maximum column abundance is 
about 5.5 x 10 9 cnr 2 ; whereas in the southern hemisphere, K is 
most abundant at high latitudes, where the column abundance is 
about 4 x 10 9 cm” 2 . Because the Caioris basin was visible in the 
northern hemisphere above 25“ latitude at the time of the observa¬ 
tions and the K was weakest there, we conclude that the Caioris ba¬ 
sin was not an enhanced source of K to the exosphere, 

22. Observations of magnesium 

Magnesium was observed in Mercury’s exosphere by the MES¬ 
SENGER MASCS instrument during the second Mercury flyby 
(M2) by its 285.3-nm emission line (McClintock et al„ 2009). MAS¬ 
CS is a scanning grating spectrometer, so that individual emission 
lines are observed sequentially (McClintock and Lankton, 2007), 
Herein we analyze observations taken during the "tail-sweep,” 
“fantail," and terminator regions, respectively, as discussed in 
McClintock et al. (2009). The tail-sweep geometry is illustrated 
in Fig. 4. Because MASCS observed across the nightside 
shadow, in the analysis we correct for the fact that, although there 
is gas in the shadow, it cannot emit there because the emission 
lines are due to resonant transitions excited by the absorption of 
solar radiation, This geometry is discussed in Appendix A. The 
tail-sweep observations were designed to sweep through a vertical 
distance bounded by +3 Mercury radii (R M ) as the spacecraft ap¬ 
proached Mercury, so the spacecraft rotation angle increased 
monotonically from ±5“ to ±40“. The data analyzed in this paper 
are confined to ±4° of the equatorial plane in order to avoid issues 
involving observations at high pointing angles with respect to the 



equatorial plane. In addition to the angle north/south with respect 
to the equatorial plane, we took into account the look angle with 
respect to they-axis, which is in the equatorial plane perpendicular 
to the Sun-Mercury line, positive dawnward. These angles vary 
throughout the observational sequence. Observations in the tail re¬ 
gion are illustrated in Fig. 5A, smoothed via interpolation between 
observed points using standard Interactive Data Language (IDL) 
routines. Because of the smoothing, the image in Fig. 5A should 
not be used in a quantitative analysis but gives visual clues only. 
The physical data used in the data extraction and model are given 
in Table 3. The noise-equivalent column is 6 x TO 7 Mg atoms cm 2 
from the surface to 16,000 km behind the planet, and 4 x 10 7 from 
16,000 km to 40,000 km behind the planet, so S/N is -10 in this re¬ 
gion. The solar spectrum is shifted blueward due to Mercury’s mo¬ 
tion toward the Sun, but the flux absorbed by air atom at rest with 
respect to the planet is on the red wing of the solar line. The g-va- 
iue in Table 3 reflects this frequency shift (Kilien et al„ 2009). 

A least-squares fit to the data was run on a series of Chamber- 
lain exosphere models solved at each observation point for the sec¬ 
ond Mercury flyby observational geometry for a series of 10 values 
of surface density (no) and 10 values of temperature to determine 
the minimum variance solution (see Appendix A), The tangent col¬ 
umn observed by MASCS within 4“ of the equatorial plane was fit 
to a Chamberlain exosphere (Chamberlain and Hunten, 1987) at 
the observed angles with respect to the planetary normal, taking 
the shadow into account. We also included a possible dawn-dusk 
asymmetry multiplicative factor in the solution to account for a 
possible dawn-dusk asymmetry, but for magnesium this factor 
converged to unity. The full Chamberlain exosphere model was 
used in this analysis, even in the tail region, and is appropriate be¬ 
cause, unlike sodium. Mg is not subject to radiation pressure and 
has a very long lifetime against photo-ionization, 57.5 h (Table 3). 
The best solution for n 0 and T is given where the least squares devi¬ 
ation from the solution for all data points is minimum. The magne¬ 
sium data taken in the M2 tail sequence (Fig. SB, asterisks) 
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Fig. S. (A) Observations of the magnesium tail of Mercury taken during the second MESSENGER flyby (Mcdintock et si., 20G9). The MASCS instrument is a scanning grating 
spectrometer, so that different emission lines are observed sequentially. The pattern observed in the image is merely the result of smoothing via interpolation between 
observed points using standard iDi, routines. The actual observed emission is peaked in the equatorial plane, but because the observations were taker, through die anti- 
sunward shadow, the observed emission includes only the illuminated portion of the exosphere outside the shadow. (B) A least squares solution was fit to a matrix of 
Chamberlain exospheres of different temperatures and near-surface densities. Hie models are for integrated emission across the tail, taking into account the lack of emission 
in the shadow. The best fit to the data, shown by asterisks, was for a temperature of 50,000 K and a density, n* of 16 cm s ; however, this fit is probably influenced by the 
single observation at a radial distance 2 fi M behind the planet. We allowed for the possibility of a different density on the dawn and dusk sides, but the solution converged to 
the same number density, giving no dawn/dusk asymmetry. A number of models are shown for temperatures of 10,000,15.000,25,00(3 and 50,000 K, respectively, each for 
surface number densities of 10,20, and 30 Mg atoms cm 1 . Note that we are seeing only the escaping component of the exosphere in th is region. The models are such that a 
higher surface number density will compensate for a lower temperature source, but lower temperature solutions require higher emission near the planet, than observed. More 
data points in this region would give greater confidence in the solution. 


Table 3 

Physical data for Mg. 



g-value at M2 for 285.296 nm line (vacuum) 0.3177 photons s ' atom' 1 
Photo-ionization rate at Earth 5.66 x 10 7 s' 1 

Photo- ionization rate at Mercury, M2 4627 x I0" 6 s' 1 

ionization lifetime at Mercury, M2 57.5 h 


converged to a Chamberlain exosphere with a near-surface density 
rto » 16 cm 3 and a temperature of 50,000 K (Fig. SB). However, the 
solution is by no means unique as shown by the curves in Fig. 5B 
for solutions with temperatures ranging from 10,000 K to 
50,000 K and surface number densities of 10, 20 and 30 cm 3 , 
respectively. A low temperature can be compensated by a high 
density at the surface to fit the data in the tail region. The best- 
fit solution is forced by the single low data point at 2 R M . The 
presence of an extended Mg exosphere cannot be due to radiation 
pressure, and loss by photo-ionization is negligible. 


We performed the same type of solution on the Mg “fantail” 
data obtained as the spacecraft rolled through 180° from boresight 
pointing dawnward, through north, toward a duskward line of 
sight as the spacecraft moved toward the planet from approxi¬ 
mately 2550 km to approximately 350 km behind the planet 
(Fig. 6B). This geometry was taken into account in the analysis. 
The Mg fantail data are illustrated by the asterisks in Fig. 6B, along 
with the least-squares solution shown by the solid red line, and 
models at temperatures of 25,000, 40,000, and 50,000 K, respec¬ 
tively, at surface densities, n 0 . of 5 and 9 cm 3 . The lower set of 
lines are those run at n 0 = 5 cm' 3 , and the upper set of curves are 
those nm at n 0 = 9 cm 3 . 

The solution for the fantail converges to a lower temperature, 
25,000 K, and a lower surface density, 9, than the solution for the 
tail data, but the temperature is poorly constrained in either case. 
However, the temperature derived from the fantail data is still 20 
times higher than that indicative of a photon-stimulated desorp¬ 
tion (PSD) source, which would give 1200 K (Madey et al„ 1998), 
and 5 times higher than that expected from a hypervelocity impact 
source, which would give about 5000 K or less (e.g„ Berezhnoy and 




Distance from Mercury Center (Mercury radii) 


Fig. 6. (A) Observations in rbe ‘'fantaii" were taken during the spacecraft roll beginning with the line of sight (LOS) toward dawn and ending with the LOS toward dusk as the 
LOS swept through north, There is a slight increase in emission to the northwest (dusk side), but the line of sight was closer to the planet on that side. (B) Fantail data 
illustrated by the asterisks are plotted with a least-squares solution (red), at a temperature of 25.000 K.and surface number density of 9 cm" 3 . Other models are shown for 
temperatures of25,000.40,000, and 50,000 K with surface number densities af 5 cnr s (Sower set of lines) and 9 cm 3 {upper set of lines). The plotted column abundances are 
observed line-of-sight abundances. This geometry does not constrain temperature as a result of the fact that the observations are in the near-tail region at a very limited range 
of distances from the planet. 
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Kiumov, 2008). The extreme temperature fit in the tail and fantail 
region is inconsistent with a PSD source for the Mg atoms reaching 
this region of the exosphere and probably requires a secondary 
process such as dissociation of a molecular source that was initially 
vaporized or sputtered from the surface, because the atoms are 
more energetic than they would be if derived from PSD, impact 
vaporization, or ion sputtering. This inference is consistent with 
that made from observations of Ca (Kiiien et a!., 2005). 

After the spacecraft completed the roil maneuver, MASCS exe¬ 
cuted Mg observations in the near-terminator region while still 
on the nightside (McClintock et al., 2009, Fig. 2). For this dataset 
the boresight was looking at the nightside surface, but only that 
portion of the exosphere above the shadow was illuminated and 
therefore seen in emission. We modeled the emission between 
the spacecraft and shadow, taking into account the boresight an¬ 
gles with respect to the surface. A fit to the Mg data in the near-ter- 
minator region (Fig. 7B) was less successful than the fits to the tail 
and fantail data, probably because the exosphere is non-uniform in 
this region and the data cannot be made to fit a single model. In the 
terminator region the observed line-of-sight column, that between 
the spacecraft and the shadow, sampled a restricted range of the 
exosphere. Each data point (asterisks) represents not only a differ¬ 
ent spacecraft distance but also a different distance behind the ter¬ 
minator. Fits are shown for a temperature of70,000 K (green) and a 
temperature of 10,000 K (red) with surface number densities of 50 
(lower curves) and 100 cm 3 (upper curves). As the spacecraft ap¬ 
proached the terminator, the data were fit by a progressively lower 
temperature and a higher surface number density. The data can be 
fit with a temperature of 10,000 K provided that the surface num¬ 
ber density is above 50 cm ' 3 , but the fits are non-unique. The first 
three points were not included in the fit because they represent a 
very short line of sight in sunlight. 


3. Data analysis 

3.1. Surface composition 

The surface composition of Mercury is not known at present. 
Table 4 shows abundance of elements from lunar soils that will 
be used to constrain possible surface compositions at Mercury: 
Luna 16 and 24, Apollo 11 low K and high K, and Apollo 17 orange 
glass (Heiken et aL 1991). In the following discussion, we consider 
the implications for exosphere models of a Mg abundance in the 
regoiith of 5 % by number (consistent with the mean of the Apollo 
11 and Luna measurements) and 8.6% by number, consistent with 
Apollo 17 orange glass. 

3.2. Sodium and potassium 

The potassium map obtained on January 17. 2008, at the 
McMath-Pierce solar telescope had 78 pixels in the illuminated 
portion of the planet. Zenith column abundances, along the norma! 
to the planet surface, were calculated from observed line-of-sight 
abundances for a Chamberlain exosphere. Although this is not an 
exact model, it does remove the line-of-sight effects given the 
assumption of a locally uniform exosphere. The average zenith col¬ 
umn abundance of Na in the northern hemisphere was 4.4 x 10 !0 , 
and that in the southern hemisphere was 3.3 x 10 10 ; the average 
zenith K abundance in the northern hemisphere is 7.5 x 10 s , 
whereas that in the southern hemisphere is 1,3 x 10 s cm 2 . Since 
the noise-equivalent column for K is about 7x10®, most of the K 
measurements in the northern hemisphere are consistent with 
noise, so only an upper limit is determined except for an area near 
the equator. From the zenith column abundance ratios per pixel, 
the average Na/K ratio in the north, excluding pixels for which K 


£ 



Fig. 7. (A) Illustration of the geometry of the near-dawn terminator observation. The arrows show the direction of the UVVS line of sight during each observation, and the 
coior represents the observed intensity. The observations were approximately in the equatorial plane. {8} Models fit to the near-terminator data indicate a best-fit solution at 
a temperature of35,000 K and density n (! of 75 cm” 3 ; however, we show models at temperatures of 10,000 K and 70,000 K, respectively, each at surface number density of 50 
and tOO cm" 3 . Each individual data point («) is obtained at a different distance behind the terminator, with those at s larger spacecraft distance from the planet closer to the 
terminator. As the spacecraft moves toward the dawn terminator a longer path length is observed, and the path intersects less shadow. The noise-etjuivaient column is 
6 x 10 s cm 3 from the surface to 16,000 km behind the planet, thus S/N is ~1.7 for the least emission. Points below ! .15 R M were not modeled because the path length in 
sunlight was too small. The best-fit temperature decreases and the best-fit surface number density increases as the spacecraft approadies the terminator (dawn). There does 
not appear to be a single solution in this region. A simple Chamberlain exosphere with n ( , - 75 cm 3 and a scale height of-3100 km (corresponding to 35.000 K) wouid have a 
column abundance of 23 x S0 !c atoms em u 


Table 4 

Model abundance of elements for Mercury surface materials and abundances in selected lunar soils. 

Species ■ Model 1 lanjti$-#as ■ Itfna 24;#SI\ , Apo8&li high K Apollo 11 lowK#% _ Apajfo 17 orange glass #% 

Ca- 2.15(1.3) - 4.9-- 4.7 4.4 4.44 3.1. 

Mg 5-46 ■' S3 5.9 42 4,0 .8.6. 

Na - 0.34 (035) 0.27: 022 0.39031 6.28 

K 0.03 (0.02) 0.06 0.18 0.03 0.045 
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is at noise-equivalent column, is 48.5. The average Na/K ratio in the 
south, also excluding pixels for which K is noise-equivalent col¬ 
umn. is 22 where the S/N in K is 5-10. 

Note that historical data records give the ratios of line-of-sight 
column abundances. To facilitate comparisons, we compute these 
ratios as well even though they have line-of-sight effects. The 
line-of-sight average Na abundance was 9.9 x 10 to in the northern 
hemisphere and 8.2 x 10 ,a in the southern hemisphere. The aver¬ 
age K line-of-sight abundance in the northern hemisphere was 
1.8 x 10 9 t and that in the southern hemisphere was 
2.6 x 10 9 cm 2 The average line of sight values are roughly twice 
those of zenith column values. The very large Na/K ratio in die 
northern hemisphere reflects the fact that the K abundance is very 
low in the northern hemisphere at this time. More interesting is 
the spatial distribution of K. In the northern hemisphere K is al¬ 
most entirely confined to within 30 s of the equator, whereas in 
the southern hemisphere K ts mostly south of 30°S. 

Berezhnoy and Klumov {2008} concluded that for an impact- 
vaporization source of the exosphere of Mercury, the ratio of Na 
atoms to all Na in the exosphere (atomic and molecular) is 0.7, 
whereas the ratio of K atoms to all K is 0.4. If we assume that these 
ratios hold for all source processes, and the average ratio is most 
representative of source processes, then the Na/K ratio in the exo¬ 
sphere can be approximated as 

Nap, Na su[f WNa) Na ws K ati , 

K gJS " K surf t iom2 (K) Na tMa i K S m 

where is the measured exospheric ratio (37), Nagar/Na,^! 

- 0,7, and Kgaa/Ktoti = 0.4 from Berezhnoy and Klumov (2008). Given 
photo-ionization lifetime of Na from Huebner et aL (1992) and that for 
K from Fulle et al. (2007) (Table 2), we derive tire ratio of abundances 
in the regolith, Na. 


37 =■ N-w/Karf/3.75 ■ 1.75) 


(2) 


This last expression yields a ratio of Na/K in the soil of approx¬ 
imately 5.6. For the ratio of Na/K found in the southern hemi¬ 
sphere, where the K abundance is above the noise-equivalent 
column, the observed atmospheric Na/K ratio is 22, yielding a 
Na/K ratio in the surface of 3.4. The K content of Apollo 17 pristine 
orange glass is 0.045% by number (625 pg/g). The Na/K ratio in 
Apollo 17 pristine orange glass is 6,2, and that in Luna 16 soil is 
4.5. The average of our derived ratios of Na/K in the soil is 4.5 
and is therefore consistent with the abundance ratio in Luna 16 
soil, for which the K abundance is higher than that in Apollo 17 or¬ 
ange glass but lower than that for high-K Apollo 11 lunar samples. 
The Na/K ratio derived from the northern hemisphere is consistent 
with Apollo 17 orange glass, a pyroclastic glassy deposit. 

Photon-stimulated desorption has been shown to be an impor¬ 
tant, if not dominant, source of alkalis in the lunar and mercurian 
exospheres (e.g., Madey et al„ 1998; Madey and Yakshinskiy, 1997; 
Yakshinskiy and Madey, 1999, 2004, 2005), Potter and Morgan 


Table S 

Observations of Na/K and relevant physical parameter's. 


Date- TAM Phase (•) 


December S.1B90- 0,365 274.5 77 

December j, 1990 0,359 2705 SD 

Decembers,. 1390 ' 0353 284.2 S4 

Decembers, 1980 0.348 289.1. 88 

December 10,1990 0342' 294.1 93 

jaraiary 11,1998 0.424. 129.5 82 

April 26,1999 0.444 216.S 77 

May 2,1999 . 01419 2353 66 

January 17,2008 0.334 300.2 63 


8 (fan s') 


+7.8 

-6.0 

-S.7 



(1997) showed that the variations of Na and K in Mercury's exo¬ 
sphere are correlated and depend on the flux of solar ultraviolet 
(liV) radiation. Their ratios of Na/K in the timeframe December 
6-10, 1990, varied from 130 to 260; these are much higher than 
found here but probably included line-of-sight effects and also 
may indicate that low values of K abundance in the exosphere 
are consistent with noise-equivalent column. Potter et al. (2002) 
found an average Na/K ratio in Mercury’s exosphere of about 
100, There were three dates reported in that paper on which the 
Na/K ratio was most similar to the one found here; January 11, 
1998 (67.7), April 26,1999 (66.3), and May 2,1999 (69.5). We have 
included the observations of Potter et al, (2002) and our observa¬ 
tions in Table 5 along with orbital distance from the Sun, true 
anomaly angle, heliocentric radial velocity, g-value of Na(D| + D 2 ), 
Na/K, radiation acceleration (cm s 2 ), and F10.7. The FI 0,7 index is 
derived from the 10.7-cm radio flux, which has been shown to be 
correlated to the solar Lyman a flux and is a standard proxy for so¬ 
lar ultraviolet (UV) flux. The average value of F10.7/(Na/K) is 14.5, 
regardless of whether FI 0.7 is high or low. The values of Na/K ratio 
vs, Lyman at flux normalized to the heliocentric radial distance 
(Fig. 8) give a correlation coefficient of 0.66, which implies a weak 
relationship between UV radiation and the Na/K ratio in the exo¬ 
sphere (Fig. 8). The values of Lyman a flux (photons cm 2 s ’) in 
Table 5 are from the laboratory for Atmospheric and Space Physics 
(LASP) Interactive Solar IRradiance Datacenter (L1SIRD) website 
(http://Iasp.colorado.edu/LISlRD/) (Snow et al., 2005) and have 
been rotated to the orbital position of Mercury for a solar rotation 
rate of 27 days (disk average) and the phase angle given in Table 5. 
We conclude that the solar UV flux provides an important source of 
alkalis to the exosphere but cannot be the only source. The depen¬ 
dence on UV flux for the Na/K ratio implies that Na and K vary dif¬ 
ferently with UV flux. 


3.3. Magnesium 

3.3.1. Impact vaporization at the time of M2 

The magnesium observations reported here were obtained dur¬ 
ing the second MESSENGER flyby on October 6, 2008. Mercury’s 
orbital distance was 0.342 AU. We estimated the impact vaporiza¬ 
tion rate on this date using the planar impact approximation meth¬ 
od (Meiosh, 1989, chapter 4; Morgan and Killen, 1998; Killen and 
Sarantos, 2004), and we scaled the impact flux from Earth orbit 
measured by Love and Brownlee (1993) to Mercury orbit following 
the mathematical formulation in Cintala (1992). We approximated 
the thermodynamic properties of the meteoroid impactors and the 
regolith using the "regolith” constants in Cintala (1992) and the 
impact parameters for Al onto enstatite given in Meiosh (1989), 
Appendix B. Using this formulation the total vaporization rate at 
Rorb» 0.342 AU is 2.7 x 10 7 cm~ 2 S" 1 for impact of aluminum (a 
proxy for asteroid composition) onto regolith, and the Mg 


Na/K Radiation acceleration (Na) 


214 1692 



130 181.3 

66.3 64.5 

69.5 110.7 

SO. 180.0 


FtCt? gg£ 4'ax 10” Ly oe/K 2 x 10” 

2277 50.6 5.74 43.1 

2288 14.8 5,78 44,8 

2304 8,9 $.81 46.6 

.2374 15.1. 5.82 48.1 

2407 18.5 5.82 49.8 

1017 15.0 3:91 21(7 

1058 15.9 4.18 21.2 

737 122 3.62 32.4 


Notes: H^a, - orbital distance in AU. TA -1 


1 anomaly angle. R - h.eiic 


ric radial velocity (km: 
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Fig. *. The Ns/K ratios from Porter et a!. (2002} are plotted vs. Lyman along 
with that from, this paper {data point at Ly s/id* - 32). The Na/K ratio is moderately 
correlated with the Lyman a radiation. We plot the average Na/K ratio over the 
illuminated disk to avoid complications regarding scale-height effects off the disk 
and added noise caused by iow-signal K observations as the K decreases away from 
the planet. Although a linear correlation coefficient of 0:66 does not represent a 
strong correlation, it indicates that 66% of the variation can be explained by a linear 
correlation between the ratio of Na/K and Lyman at flux. However, our data point is 
much lower than those found earlier because we rejected K values that have SjN < i, 
This indicates that although OV radiation is involved, other factors are effective in 
the ejection and loss processes. 


vaporization rate is 2.3 x 10 6 atoms cm' 2 s~ l for an 8.6% surface 
abundance of Mg or 1.3 x 10® cm 2 s"' for a 5% abundance of Mg. 

In addition to the Cintaia (1992} and Morgan and Killen (1998) 
analyses of meteoroid impact rates onto the surface of Mercury, 
several new estimates have been published. Notably, Borin et al. 
(2009) have estimated a meteoroid impact flux onto Mercury 
170 times the Cintaia rate. This of course has important implica¬ 
tions concerning the exosphere. Killen et al. (2001) used the Cintaia 
(1992) mathematical formulation for transfer of material from 
Earth orbit to Mercury orbit with the Love and Brownlee (1993) 
interplanetary dust particle (1DP) flux rates and the Morgan and 
Killen (1998) impact vaporization code to estimate the magnitude 
of various contributions to the sodium exosphere and obtained an 
estimate of 20-30% of the sodium exosphere contributed by im¬ 
pact vaporization. Thus, an increase of the Love and Brownlee 
(1993) flux by a factor of only 3-5 would supply the entire exo¬ 
sphere. Morgan et al. (1988) estimated that meteoroids supply 
up to 1.4 x 10 24 Na atoms s 1 to the exosphere, or a factor of about 
5 greater than the Cintaia rate, and concluded that meteoroid 
vaporization can supply the full Na exosphere of Mercury. The 
main reason for the larger rate given by Morgan et al. (1988) is 
the use of a velocity distribution function from Southworth and 
Sekinina (1973). An impact flux 170 times the Cintaia rate would 
either overwhelm the exosphere, or we would have to conclude 
that impact vaporization is either less efficient than previously 
estimated, that impact vapor quickly condenses, leaving only 2% 
of the vapor in the exosphere, or that the vaporization products 
are not fully measured, possibly being in molecular form. The dis¬ 
tribution of various elements in the exosphere appears to be incon¬ 
sistent with uniform impact vaporization by micrometeoroids (e.g., 
Burger et al, 2010), mainly because of high-latitude enhance¬ 
ments, a dawn enhancement in Ca, and rapid variability. 

The Mg observations in Mercury’s tail region were fit by a 
15.000-50,000-K gas, which is much hotter than would be ex¬ 
pected from direct impact vaporization. Although the temperature 
is not well constrained, and indeed the velocity distribution is 
likely to be non-thermal. it is dear that a hot vapor is required to 
fit the data. Very hot (30,000 K) Ca atoms have been observed in 
Mercury’s exosphere from ground-based observations (Bida et ai„ 
2000; Killen et a!., 2005) and were attributed to vaporization of cal¬ 
cium in molecular form with subsequent dissociation, which im¬ 
parts additional energy to the atomic fragments. We suggest that 


the same process produces the hot Mg that we have observed in 
the tail region. The dissociation energy of MgO is 3.35 eV (Murthy 
and Bagare, 1978). If dissociated by a Lyman a photon (10.2 eV) the 
products would have an excess energy of 6.8 eV between them, if 
each fragment lias 3.4 eV, then the equivalent temperature is 
40,000 K. If the gas temperature is 40,000 K, the mean velocity of 
the gas is 5.2 km s which is above escape velocity (4.2 km s" 1 ). 
Thus those Mg atoms that are derived from dissociation of MgO 
will escape. Their lifetimes will be the photo-ionization lifetime 
(Table 6). The Mg ionization rate is 5.66 x 10 7 s”' at Earth (quiet 
Sun) or 4.84 x 10 6 s ' 1 at Mercury during the second MESSENGER 
flyby. Thus the Mg photo-ionization lifetime was 2.0 x 10 s s, or 
57 h. At a velocity of 5.8 km s' 1 , a Mg atom would travel 475 plan¬ 
etary radii on average before becoming ionized. The Hiil sphere, 
where the gravitational field of the Sun dominates over that of 
Mercury, is at ~81 R M - Therefore ionization is not a major loss pro¬ 
cess for Mg in the region that we observe. Note that even though 
the Mg atoms in the tail have very high velocities relative to the 
planet, the g-values. and hence the derived column abundances, 
only differ by 5% between that for a Mg atom at rest with respect 
to Mercury (g-value* 0.317) and one traveling at 5 km s 1 anti- 
sunward (g-value=0.300). Therefore the assumed g-values do 
not introduce a large error in the derived column abundances even 
if the assumed velocities are incorrect. 

The fraction of Mg(atoms)/Mg(total) from impact vaporization 
at quenching was estimated by Berezhnoy and Klumov (2008) to 
be less than that of the alkali species. At 4000 K. the equilibrium 
fraction of Mg to MgO is about 2 (Berezhnoy and Klumov, 2008). 
We assume that half of the Mg + MgO will condense, leaving a 
gas release rate of 0.7 x 10 s < dnjdt < 1.2 x 10 6 for Mg abundances 
of 5% and 8.6%, respectively. If we assume that 30% of the impact- 
derived Mg is in the form of MgO, then the very hot portion of the 
Mg exosphere should have an impact vapor source rate of 
0.3 • 1.2 x 10 s on' 2 s ! » 3.6 x 10 s cm -2 s 1 for an 8.6% Mg abun¬ 
dance in the regoiith or 2 x 10 s cm 2 s 1 for a 5% abundance, using 
the vapor code described in Morgan and Killen (1998). The ob¬ 
served column is estimated to be 7 x 10 9 cm 2 norma! to the sur¬ 
face, If the exosphere is that portion from the surface to the Hill 
sphere (81 planetary radii) and the Mg atoms have a mean velocity 
of 5 km s \ then the lifetime of the Mg atom in the exosphere is 
4 x TO 4 s. Then if the production rate of Mg is 6.75 x 10 5 s ’, and 
the lifetime is 4 x 10 4 s, we would expect to see a column of 
2.7 x 10 10 Mg atoms. Given that we see 26% of that, we conclude 
that this is the fraction of the vapor initially in the form of MgO 
that was dissociated to produce the very hot atomic Mg observed, 
This inference is consistent with the result from Berezhnoy and 
Klumov. namely that 30% of the vapor is in the form of MgO. We 
conclude that impact vaporization of Mg-bearing minerals can pro¬ 
duce the observed Mg tail given an abundance of Mg in the regoiith 
of 5-8.6%, if roughly half of the vapor condenses and 30% of the va¬ 
por is in the form of MgO as estimated by Berezhnoy and Klumov 
(2008). This relies on the assumption that the observable lifetime 
of Mg is the lifetime of the atom to the Hill sphere. 




ionization rate, s ' Lifetime at 0.342 
(at 1 AU. active Sttn-j AU (s) (quiet Sun*) 


Na 5.99x10'* 6.52x10 

Mg 5.66 x 10' 7 1.08x10 * 

C4 7 JO x !G~ S 7$ x tO - 

tC 2.2 x 10 5 24 x 10 5 


1.95 x 30 4 
2.07 x UP 
1,67 x 10- 
53x10? 


J ionization rates depend on the solar i 
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3 32. Ion sputtering at the time of M2 

Sputter yields of Mg from MgO by low-energy He and Ar ion 
beams were measured as a function of ion energy by Hine et aL 
(2007) and fit to the following equations: 

Y H '^a(y l/Eo-l) (3) 

where V H , is the MgO sputter yield from incoming He ions, 
a =0.191, and E 0 “241 eV; 

Y* - 0(1 - v%]tf a (l + *(v^rn)) ( 4 ) 

where y Ar is the sputter yield of MgO by impacting Ar ions, 
a = 0.0644, b-7.61, and E 0 « 36.5 eV. These measurements were 
made for E < 460 eV; however, we will employ the equation to esti¬ 
mate yields by solar wind ions. The sputter yield of Mg by 3.3-keV 
He ions in the solar wind would therefore be about 0.5/ion. 

Ignoring for the moment highly charged ions, which normally 
comprise less than T % of the solar wind, and assuming that He ions 
are 4% of the solar wind, the solar wind yield would be 

^1 = 0.02/^ (5) 

where/ sw is the solar wind flux onto the surface of Mercury and fug 
is the number fraction of Mg in the soil. We estimate that the solar 
wind flux onto the surface of Mercury during the second flyby was 
10 25 ions s 1 (e.g„ Benna et al., 2010), although this may be an opti¬ 
mistic estimate (e.g„ Sarantos et al., 2007). If we assume that the 
fraction of Mg in Mercury’s regolith is 8%, then the sputter yield is 
about 1.6 x 10 Z2 atoms s" 1 or a planet-wide average of 
2 x 10 4 Mg cm 2 s ' 1 . If the lifetime of the Mg atom to the Hill 
sphere is 4 x 10 4 s, the column from sputtering would be about 
8 x 10 8 Mg cm 2 , or about 10% of the observed column. Although 
the sputter yield averaged over the planetary surface is two orders 
of magnitude less than the estimated hypervelocity impact yield, 
the local yield in areas of open fields may be comparable, especially 
at times of enhanced solar wind flux. The local proton flux to the 
surface can be as high at 10” ions cm" 2 s 1 and averages 10 8 in 
the open areas (Benna et al„ 2010). With a yield per ion of 0.02 
/Mg or ~0,002 f svn the local sputtered Mg rate averages about 
2 x 10 s cm 2 s 1 in the cusp regions during the first MESSENGER 
flyby when the magnetosphere was quiet and the interplanetary 
magnetic field (IMF) was northward, impeding solar wind access 
to the surface. The sputter yield may be as high as 0.2/ion, in which 
case the sputtered Mg rate would be about (1 -2) x 10® cm 2 s" 1 in 
the cusp regions during the first MESSENGER flyby. This local rate 
would be comparable to the impact vaporization rate and would 
produce hot atomic Mg. If this ion-sputter yield adds to the impact 
vaporization yield, then the local yield could produce local enhance¬ 
ments equal to 10 times the observed average column, as observed 
by MESSENGER during the third flyby (McClintock et al., 2009). The 
estimated sputtering source is therefore a small fraction of the total 
Mg exosphere but could be important for populating the tail region. 

Because the disturbed solar wind can be comprised of enhanced 
fractions of heavy elements, for instance up to 30% He ions (Wang, 
2008), sputtering by highly charged heavy ions may play an impor¬ 
tant role at Mercury (Kallio et ai.. 2008). Sputter yields by highly 
charged ions are known to be significantly enhanced over those 
by singly charged ions, especially in insulators. This effect is known 
as "potential sputtering" and has been extensively discussed in the 
literature (e.g., Kallio et ai., 2008). Although the charge-state effect 
is not as great in highly ionic oxides such as MgO, there is still a sig¬ 
nificant charge-state effect (Hayderer et al., 2001). If we can corre¬ 
late enhanced ion sputtering with solar energetic particle (SEP) 
events at Mercury, we may be able to determine whether the 
highly charged ions are an important source of sputtered atoms. 


Additional laboratory work on sputter yields by highly charged 
ions is needed to more fully understand this process. 

3.3.3. Estimate of photon-stimulated desorption of Mg 

Magnesium can be desorbed from MgO by 265-nm laser pulses 
(Joly, 2002), and PSD yields of Mg from MgO are comparable to 
those for Na. Although MgO comprises about 15-26% of pyroxene, 
and the pyroxene abundance is about 40-60% for mare basalts on 
the Moon, the Mg is present as a silicate. Therefore the PSD yields 
for Mg from Mercury's surface minerals is highly uncertain at the 
present. Mg would be desorbed by PSD only if a significant amount 
of MgO were re-deposited on the surface as a result of impact 
vaporization. 


4. Discussion 


Potter et al. (2002) found that the ratio of sodium to potassium 
in the exosphere of Mercury was highly variable and averaged 
100 ± 20, much larger than values observed elsewhere in the Solar 
System (Table 7). Hie Na/K ratio in the lunar exosphere is close to 
that in the lunar crust We will argue that the Na/K ratio in Mer¬ 
cury's exosphere is in fact much larger than that ratio in the rego¬ 
lith, and that very large ratios of Na/K in Mercury’s exosphere may 
be due to very low abundances of K in the exosphere, where mea¬ 
surements are dominated by noise. 

The Na/K ratio in the exosphere of the Moon is in the range 4.4- 
6, not far from overall average ratio of Na/K in the surface lunar 
rocks, which is in the range 7-9. Potter et al, (2002) found that 
radiation acceleration losses were not sufficiently different for so¬ 
dium and potassium to explain the Na/K high ratio. They proposed 
that the photo-ionization loss and recycling rates were sufficiently 
different for sodium and potassium to account for the high ratio. 
The higher mass of the potassium atom relative to sodium (39 
vs. 23) results in a smaller scale height and a larger gyroradius 
for potassium relative to sodium, which may result in more rapid 
net loss of potassium from the planet. However, Sarantos (2005) 
considered relative loss rates of Na and K photoions to the solar 
wind and concluded that the difference was too small to explain 
die variations in the Na/K ratio in the exosphere. Our results show 
that the Na/K ratio in Mercury’s exosphere is in the range 22-49, 
where the ratios are taken only where the observations are above 
the noise-equivalent column. 

A possible explanation for the large variation in the published 
Na/K ratios at Mercury is the use of improper g-values. This is cer¬ 
tainly the case for the high value of Na/K reported in the Earth’s 
atmosphere (Gault and Rundle, 1969). Although the g-vatues were 
scaled for the radial velocity of Mercury at the time of the observa¬ 
tion, there are two effects that were not taken into account. First, 
the radial velocity not only affects the continuum at the rest fre¬ 
quency of the atom in the atmosphere of Mercury, but it affects 
the multiplet structure, since each transition is independently af- 
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fected by the structure of the solar spectrum. Secondly, the velocity 
of the atom inside the atmosphere was not taken into account in 
the data reductioa The radial velocity of the bulk atmosphere 
may vary with the solar zenith angle and the viewing geometry. 
Third, the K emission is weak and the uncertainty in the abundance 
is large. Nevertheless there are valid explanations fora variable N a/ 
K ratio in the exosphere of Mercury. The questions here are why is 
the source rate variable and why is it different for Na and for K, 
There is a moderate correlation between Na/K and Lyman at 
flux, with a correlation coefficient of 0.6. The measured ratio is 
either large (on the order of 150-250) for Lyman at greater than 
5 x 10' 1 photons cm 2 s or small {on the order of <60} for Lyman 
a less than 5 x 10" photons cm 2 s"'. One explanation is that the 
rates of supply of volatiles to the extreme surface, or the upper¬ 
most few monolayers, and thus to the exosphere, are diffusion-lim¬ 
ited, and the diffusion rate for K is less than the diffusion rate for 
Na (Killen and Sarantos, 2004). Diffusion of atoms from the interior 
of grains to the extreme surface, or the topmost few monolayers, is 
required in order to repopulate the layers from which PSD and ion 
sputtering can eject atoms (Killen and Sarantos, 2004). Diffusion- 
limited supply is a well-known phenomenon in sputtering (e.g., 
Roth, 1983). The rate at which Na and K diffuse to the surface of 
grains of sizes 1 pm to 1 m, at temperatures found on the surface 
of Mercury, was calculated by Killen and Sarantos (2004) using dif¬ 
fusion data from Nyquist et al. (1979). Note that this is a funda¬ 
mentally different problem than that considered by Cheng et al, 
(1987) and by Killen and Morgan (1993,1994), who treated Knud- 
sen flow from depth, Killen and Sarantos (2004) showed that the 
volatiles initially near the rims of grains are rapidly depleted, and 
that the rate at which a trace volatile is delivered to the extreme 
surface of a grain, where it can be acted upon by either photons 
or ions, is not only a function of temperature but also of grain size 
and time since the grain was first exposed to space weathering, (n 
other words, there is a maximum ejection rate for a particular ele¬ 
ment set by the rate at which that element can diffuse to the ex¬ 
treme surface of the grain where it can desorb, and this rate is a 
function of grain size, phase (crystalline or amorphous), tempera¬ 


ture, and weathering. We call this the diffusion-limited supply 
rate. Of course, meteoroid impacts can vaporize entire grains or 
boulders, depending on the impactor size, and this must be fac¬ 
tored into a global solution. 

Fig. 9 shows that for a diffusion coefficient D » 10" 15 (cm 2 s ') 
all grains smaller than about 0.1 cm in radius will be able to supply 
the atmosphere at the required rate given by the known exosphere 
loss rate, assumed here to be 10 7 cm 2 s', until half of their solute 
is lost For progressively larger grains, a lesser fraction of solute 
will be lost at the required rate, and thereafter the rate will drop 
exponentially. This result is completely consistent with that of Kil¬ 
len (1989), who determined that diffusion out of glass grains is too 
fast and diffusion from a semi-infinite slab of crystalline mineral is 
too slow to maintain the exosphere. This is because small, hot 
grains, and glass grains, lose their volatiles rapidly, and large, crys ¬ 
talline grains or a slab lose their volatiles at a much slower rate. For 
smaller values of the diffusion coefficient, progressively smaller 
grains will be unable to deliver their solute to the extreme surface 
at the required rate. One-centimeter-radius spheres can maintain a 
maximal thermal vaporization rate at the subsolar point (at perihe¬ 
lion) if they are glass. At aphelion, even mineral grains can main¬ 
tain the maximal thermal loss rate for the lifetime of the grain 
on the surface. 

Given a fractal distribution of grain sizes, arid a distribution of 
temperatures on the surface, on average the exosphere is main¬ 
tained at the observed rate (Killen and Sarantos, 2004). How does 
this explain the observed variability? 

Boundaries on the diffusion coefficients for Na and K. respec¬ 
tively, are shown in Fig. 9 for glass from the leftmost boundary to 
the right-pointing arrow and for minerals from the rightmost 
boundary to the left-pointing arrow. Note that the expected K diffu¬ 
sion coefficients are in the region where the fraction of solute that 
will degas at a rate of 10 7 cm 3 s 1 is a steep function of the diffu¬ 
sion coefficient, D (and thus of temperature). Suppose that an event 
occurs that causes an increased yield of alkalis to the exosphere, 
such as a coronal mass ejection (e.g., Kallio and janhunen, 2003). 
Fig. 9 shows that Na atoms will be able to fill holes created by 
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desorption much more rapidly than will K atoms. All grains larger 
than 1 pm can supply Na at a rate of 10 7 cm 2 s \ but only grains 
larger than 100 pm can supply K at that rate. The fraction of solute 
lost per unit time at a given temperature is much greater for Na 
(pink boundary) than for K (blue boundary). The region of interest 
for Mercury is between the bold blue lines on the log diffusion coef¬ 
ficient scale for K and between the bold pink lines for Na. Thus, the 
Na/K ratio in the atmosphere may increase after an energetic event 
in which most surface atoms are desorbed because the Na sites will 
be rapidly refilled, but the K sites will take more time to recover. 


5. Conclusions 

The existence of very hot refractory species in Mercury’s exo¬ 
sphere is consistent with ejection of molecules from the surface, 
with subsequent dissociation into hot fragments. The hot Mg tail 
is consistent with a surface abundance of 5-8% Mg, if 30% of im¬ 
pact-vaporized Mg remains as MgO at quenching and half of the 
impact vapor condenses. All of the ejected MgO molecules are as¬ 
sumed to be dissociated to populate the tail region with the ob¬ 
served hot atomic Mg. Globally, ion sputtering is not a major 
source of Mg, but locally the sputtered source may be comparable 
to or greater than the impact vapor source. Photons can desorb Mg 
from MgO, so PSD may be an important source for Mg on the day- 
side if MgO is deposited on the surface by some process such as im¬ 
pact vaporization of a regolith mineral. Mg has not been observed 
on the dayside of Mercury to date, and because of the sparsity of 
data we cannot determine the global distribution of Mg at this 
time, Nonetheless, the terminator data show that the temperature 
of the Mg gas decreases and the Mg number density at the surface 
increases as the dawn terminator is approached; 

We conclude that the Na and K in Mercury’s exosphere can be 
derived from a regolith composition similar to that Luna 16 lunar 
soil in which the abundances are 0.0027 for Na and 0.0006 for K 
(Heiken et al„ 1991). The variability of Na and K in the exosphere 
is shown to be weakly correlated with the solar UV flux at Mercury 
with a correlation coefficient of 0.66. This variability of the Na/K ra¬ 
tio may be related to the differences in diffusion from the interiors 
of regolith grains to the extreme surface where they are ejected by 
PSD or ion sputtering. The variability in the exosphere is larger 
than local differences in the Na/K abundance, if they are similar 
to differences between Apollo orange glass and Luna 16 samples. 
Since Caioris basin was visible in the northern hemisphere be¬ 
tween 23° and 58 c from the sub-Earth longitude on the illuminated 
portion of the disk, and exospheric potassium was found to be 
weak in that region, we conclude that K is not enhanced over the 
Caioris basin at that time, in contrast to the findings of Sprague 
et al. (1990): Very large values of Na/K are indicative of small abun¬ 
dances of K and are influenced by noise. 

Observations of Mg in the near-terminator region appear to 
show that the surface number density increases and the tempera¬ 
ture of the gas decreases as the dawn hemisphere is approached. A 
single mode! does not match all of the data, indicating the exis¬ 
tence of a non-uniform exosphere. New models of the Mercury 
magnetosphere derived from MESSENGER observations should 
make it possible to test loss and recycling rates. The newly discov¬ 
ered element in Mercury’s exosphere. Mg, dearly is promoted into 
the exosphere by processes that are different from those ejecting 
the more volatile species, Na and K. Whereas the temperature 
of the volatile spedes, Na and K, has been shown to be approxi¬ 
mately that expected for a PSD source, or 1200 K, the temperature 
found for Mg is on the order of tens of thousands of degrees. Na and 
K exist in an extended anti-sunward tail because they are acceler¬ 
ated by radiation pressure due to the strong resonance lines in the 
visible region of the solar spectrum. There is no such strong reso¬ 


nance line for Mg in the visible region where the solar flux is 
strong; hence radiation pressure is negligible for Mg, Nevertheless 
we measured Mg in the anti-sunward tail. The conclusion is that 
the source process promoting Mg to the exosphere produces the 
hot Mg. This is very similar to the conclusion that we have found 
for Ca (Kilien et al., 2005). These refractory species are ejected into 
the exosphere by a process that imparts ~2 eV to the Mg atom and 
~3 eV to the Ca atom. We do not yet have enough data to deter¬ 
mine the global distribution of Mg. However, the data shown here¬ 
in indicate that the Mg is non-uniform. 
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Appendix A 

The comparison of data to the Chamberlain model is achieved 
by first generating model predictions of the exosphere over a range 
of temperatures and densities along the line of sight for the obser¬ 
vations. Then the model-predicted column abundances are com¬ 
pared to the data in a least squares analysis to retrieve the 
temperature and density. The analysis can be broken down into 
categories for each observation regime. 

A.1. Tail-sweep observations 

For the tail-sweep observations the spacecraft is not in the 
shadow of the planet. The line of sight oscillates up and down 
relative to the planet's equatorial plane, causing it to intersect 
the planet’s shadow for some of the observations and not 
intersect the shadow for others. A subset of observations was 
selected where the lines of sight were pointing near the equato¬ 
rial plane and passed through the shadow of the planet. Refer¬ 
ring to Fig. 4. the column abundance along the portion of the 
line of sight in the dawn-side hemisphere. S 2 , was calculated 
over a range of temperatures and densities using R 2 and <h as in¬ 
puts to the Chamberlain model and with 5-> extending from R-. to 
infinity. The column abundance along the portion of the line of 
sight in the dusk-side hemisphere, S,, was determined over a 
range of temperatures and densities by calculating the column 
abundance from R, to infinity in the direction of the spacecraft 
and from R sc to infinity and then taking the difference between 
the two. A least squares analysis was performed over tempera¬ 
ture and density for both the dusk-side and dawn-side Sines of 
sight. For a single temperature, the Chamberlain-predicted 
column abundances from the dawn-side hemisphere for the en¬ 
tire range of densities were added one by one to each of the col¬ 
umn abundances from the dusk-side hemisphere over their 
range of densities for the same temperature. This was done for 
each of the subset of observations and compared to the corre¬ 
sponding set of measured column abundances. The difference is 
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D(T, p i3m , p iu!k ) = ]T [Nj - M(T, Pima , p dusk )f (At) 

where N is the measured column abundance, M is the model-pre¬ 
dicted column abundance at a specific temperature (T), dawn-side 
density (pdaw»)> and dusk-side density (pdmk), and the sum is over 
the number of observations (n). The result is a three-dimensional 
matrix of D where the minimum value is interpreted as the best 
fit of the data to the model. 


A 2. Fantail observations 


When the spacecraft is in shadow, the line of sight can only exit 
the shadow. The column abundance was determined over a range 
of temperatures and densities along the line of sight from the point 
where it exits the planet’s shadow to infinity. The fantail observa¬ 
tions were analyzed by comparing ail of the fantail data with the 
corresponding Chamberlain-predicted column abundances over 
temperature and density. The difference is given by 

D(T, p) = [N, - M(T, p)f (A2) 


where N is the measured column abundance, Af is the model-pre¬ 
dicted column abundance at a specific temperature (Tj and density 
(p), and the sum is over the number of observations (n). The result 
is a two-dimensional matrix of D where the minimum value is 
interpreted as the best fit of the data to the model. 


A3, Tenninator observations 

When the spacecraft is not in the planet's shadow but tire field 
of view intersects the shadow and then intersects the planet, the 
Chamberlain column abundance over a range of temperatures 
and densities was calculated from the point of line-of-sight inter¬ 
section with the shadow to infinity (this column will pass through 
the spacecraft) and from the spacecraft to infinity (this will be in 
the direction opposite to the line of sight). The second quantity 
was subtracted from the first, and the same least squares analysis 
was performed as for the fantail observations. 
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